Newtonian limit of conformal gravity 
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We study the weak-field limit of the static spherically symmetric solu- 
tion of the locally conformally invariant theory advocated in the recent past 
by Mannheim and Kazanas as an alternative to Einstein's General Relativ- 
Q\ ' ity. In contrast with the previous works, we consider the physically relevant 

case where the scalar field that breaks conformal symmetry and generates 
fermion masses is nonzero. In the physical gauge, in which this scalar field 
is constant in space-time, the solution reproduces the weak-field limit of the 
Schwarzschild-(anti) De Sitter solution modified by an additional term that, 
depending on the sign of the Weyl term in the action, is either oscillatory 
CN ■ or exponential as a function of the radial distance. Such behavior reflects 

^ ■ the presence of, correspondingly, either a tachion or a massive ghost in the 

spectrum, which is a serious drawback of the theory under discussion. 

^J" ■ PACS number (s): 04.50.+h 

o : 

ON . 

In a series of papers (see and references therein), Mannheim and Kazanas explored the 
possibility that the gravity is described by the conformally invariant theory with the key 
^! ingredient in the action being the Weyl term 

S: J w = -afd 4 Xy/^gC Xlu/K C x ^ 

! = -2a f d A x^g (R^R^ - \R 2 ) + boundary terms, ^ ' 



where C X ^ VK is the conformal Weyl tensor and a is a purely dimensionless gravitational 



coupling constant (we use the system of units in which % = 1 and c = 1). In particular, they 
obtained the complete conformally static spherically symmetric solution of the theory 
Eq. ([!]) with the line element given by 

ds 2 = C 2 {x)ds 2 , dsl = -G{r)dt 2 + dr 2 /G{r) + r 2 dQ, (2) 

where C(x) is an arbitrary nonzero function of the spacetime coordinates x, and G(r) is 
given by 

G(r) = 1-/3(2- 3/? 7 ) /r - 3/3 7 + 7 r - nr 2 . (3) 

Here, (3, 7, and k are integration constants. Having tacitly assumed that test bodies move 
along the geodesies of the metric with the line element ds^ of equation (Q), Mannheim and 
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Kazanas then claimed to recover the Newtonian term (oc 1/r) in the potential of solution 
(§) of the conformal gravity theory and also suggested fl],[|,[7| that the additional linear term 
7r in Eq. (^) might account for the flat galactic rotation curves without having to invoke 
dark matter. 

However, solution (fj), (|3]) of the purely gravitational conformal theory defined by Eq. (P 
is not quite relevant to the observations, since it is obtained without regard of the matter 
part of the theory that includes the mass generation mechanism for the elementary particles 
and thereby for test bodies such as stars and planets. Such a feature of this solution is 
reflected in the unrestricted freedom of choosing the conformal factor C (x) in Eq. ([|) which 
clearly affects the timelike geodesies of the metric, but which is totally undetermined thus 
far. Moreover, the electrovac generalization of solution (|2]), (|3|) was previously obtained by 
Riegert ||, who also asserted that one of the integration constants can be eliminated by 
further coordinate and conformal transformations. This property of the solution, with 7 



being such constant, was noted also in [10,11 1 and very recently explicitly demonstrated in 



All this makes very problematic the use of the metric given by the second expression 
in Eq. (0) and by Eq. as an observable one. 

In the present paper, we consider this problem taking the matter to be represented by 
the following generic conformally invariant action @,f|,[3) 



• M 



/ d A xy/=g~ [d^Sd^S/2 + AS 4 - S 2 R/12 + i^(x)V^ - (S#] , (4) 



where ip is the fermion field, S is the scalar field, R is the curvature scalar of the metric, 
and A and ( are dimensionless coupling constants. In the theory defined by Eqs. (P, 
once the scalar field S is everywhere nonzero it can be gauged to an identical constant So 
by a conformal transformation. In this gauge, the fermion part of the action acquires the 
standard form with constant mass, hence all physical effects receive standard description; in 
particular, massive particles and test bodies move along the timelike geodesies of the metric 
as in General Relativity. It is clear that since conformal symmetry is broken and there are 
massive particles in the real world, one should take solutions with S being nonzero. The 
physical vacuum is then regarded as the state without excitations of the rest of the matter 
fields, in our case, the field ip. 

We consider solutions outside a compact source formed by the matter fields (represented 
in our model by the single field ip). The equations of the theory, first introduced by Bach 
fi"3| for the case of generic matter, have the form 

AaW^ = TV, (5) 

where the two sides stem, respectively, from the variation of the action (|I|) and the action 
(d) with respect to the metric, and the expression of the stress-energy tensor T^ u in the 
gauge S = So and with the ip field being zero is given by 



T, u = -Si (R pv - g^R/2) /6 - \S* 9/a ,. (6) 

Equations ((^) with the right-hand side given by @ are nothing but the Bach-Einstein 
equations with the cosmological constant term — the last term in Eq. @ (see, e.g., |[L4|| ). 
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Note that the left-hand side of Eq. (Q) is identically traceless, and it is convenient to rewrite 
system (H) as 

4aWV = V, R = 2AXS 2 , (7) 

where T^ u = —Sq(R^ u — g^R/A) /6 is the traceless part of the stress-energy tensor T^ u , 
and the second equation of system (|7|) is the trace of equation (||). 

We restrict ourselves to the static spherically symmetric case. As we explained above, 
we are interested in the situation where T^ u is given by Eq. (^]) with constant nonzero S . 
In this gauge, a static spherically symmetric metric can be put in the form 

ds 2 = -B(r)dt 2 + A(r)dr 2 + r 2 dtt. (8) 

Performing transformation of the radial coordinate as in 0, it is convenient to rewrite this 
metric in the form 

ds 2 = C 2 {p) [-D{p)dt 2 + dp 2 /D{p) + p 2 dVt] , (9) 

where r 2 (p) = p 2 C 2 (p), and the functions C(p) and D(p) are simply related to A(r) and 
B(r) (see 0). 

Due to spherical symmetry and stationarity, the first set of equations (|7|) will give two 
independent equations which are conveniently chosen to be 

4a (W ° - Wl) = T ° - 71 , 4aW n = %i , (10) 

where index "0" labels the time t and index "1" labels the radial coordinate p. The expres- 
sions for the quantities in the left-hand sides of equations ( |T0"D for the metric of Eq. (||) were 
obtained by Mannheim and Kazanas [l|JJ: 

1 {D'D'V D" 2 DDW-D'D" DD" + D' 2 2DD' D 2 1\ , , 
W « = 3CH){— 4 ~p — + — -y + 7')' <12) 

where the superscript "(rt)" indicates the nth derivative with respect to the radial coordinate. 
Note that the function C(p) has been factored out in equations (|TTj) and ( |T2"D due to the 



conformal symmetry of the theory. The expressions for the right-hand sides of equations 
(0) are 

12 \2D p 2 D Dp J 
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where F = logC. The scalar curvature is given by 

Q{p 2 DC')'-C{p 2 {l-D))" 

In a general case, it appears to be difficult to obtain an exact solution for C(p) and D(p). 
However, it is possible to obtain the weak-field limit of the solution. Let the physical metric 
of Eq. (H) in the spatial region of interest be sufficiently close to the flat one, so that 

A(r) = l + ea(r), B(r) = 1 - eb(r), (16) 

where e is an auxiliary small parameter to be set equal to unity in the end. Then, in the 
same region, the functions C(p) and D(p) of the metric of Eq. (§) and the radial coordinate 
transformation r(p) to the first order in e are given by 

C(p) = l + ef(p), D(p) = l-eh(p), r(p)=p[l + ef(p)}, (17) 

with the functional relation 

a(r) = h{r)-2rf'(r), b(r) = h(r) - 2/(r). (18) 



To obtain the system of equations for the functions h(p) and /(/?), we linearize the 
equations of system (|7|) for the metric of Eqs. (^), ([H]) in the small parameter e. First, we 
note that the scalar curvature R of this metric is of order e. Hence, the second equation of 
system (|^) implies that the dimensionless value of XS^r 2 should also be at least of order e 
in the spatial region under consideration. On the observational grounds, this restriction on 
the value of XS^r 2 is quite natural since this value represents the effect of the cosmological 
constant which is believed to be small on the galactic and stellar spatial scales. However, 
from the theoretical viewpoint, such a restriction constitutes the fine-tuning problem of the 
cosmological constant. The solution of this long-standing problem is absent, so we formally 
replace A by eA, thus taking into account the smallness of the corresponding parameter. 

It remains to linearize equations (|i~0|) and the last equation of system (0) in the small pa- 
rameter e. Omitting the simple but cumbersome calculations, we present the corresponding 
result in the form of the following system for the functions h(p) and /(/?): 

- {ph) {i) = 6ppf", (19) 

1 /, ^ h" 2h' 2h\ p fh" h ,.. 2/'\ , , 

U(3) + + — =^ - 3 f + -U , 20 

3p V P P P 3 J 2 V 2 p 2 p ) 

6(p 2 /')'=(A)" + 24gp 2 , (21) 

where we made the notation 
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V 



o2 

24a 



q = XS 2 . 



(22) 



We proceed to the solution of these equations. First, it is convenient to set 

h(p) = 2m/ p - 2qp 2 + v(p) , (23) 

where v(p) is the new unknown function and m is a constant. Equations (|i9|) and ([H]) can 
be integrated once. The integration constant that appears after the integration of Eq. ([19]) 
can be set to any value by rescaling the time and length which shifts f(p) by a constant. We 
use this property to eliminate the parameter q from the equations for v(p) and f(p). The 
integration constant that appears after the integration of Eq. (|2lD can be eliminated by the 
redefinition of the constant m in Eq. fl2"3]). As a result, we obtain the following system of 
equations: 



pv 



(3) + 3v" + 6p (pf - f) = 



(24) 



v 



3p 



w (3) H pv" 

p 4 



^ + (4 + W) t; + gP pr _3 pr = 



2p 3 



u' + 2v/p - 6f = . 



(25) 
(26) 



Now it is convenient to proceed to the new unknown functions a(p) and b(p) related to 
v{p) and f(p) as follows [cf. Eq. (0)]: 

(27) 



a(p) = - 2pf(p) , 6(p) = v(p) - 2f(p) . 



From system (g4| 
6(p): 



one easily obtains the following system of equations for a(p) and 
b" + 2b'/p + pb = 0, p6' + 2S = 0. (28) 



Its solution is straightforward and depends on the sign of the constant p that coincides with 
the sign of a [see Eq. (E3)]. First, we consider the case where p > 0. We obtain 



a(p) = n 



sin {kp + 0) 
P 



k cos (A;p + 0) 



(29) 



Kp) 



2nsin (kp + <p) 
P 



(30) 



where = y/p, and n and (f) are integration constants. Combining Eqs. 
(H|), (p0|), we eventually obtain the following solution of our problem: 



, p3h, (0), and 
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a(r) 



2m 



— 2qr + n 



sin (kr + <p) 



— k cos (kr + (p) 



(31) 



2m + 2nsin (kr + 0) 2 
o(r) = 2gr . 



(32) 



We see that in the Newtonian limit, apart from the universal term qr 2 , there arises the 
additional gravitational potential 



V(r) 



m + n sin (kr + 



(33) 



in which the constants m, n, and <p are to be related to the source. The constants k = ^/p 
and q are universal and are given by Eq. (|2"2"|). 

We note that the linearized static spherically symmetric solutions in a generic (not con- 
formally invariant) second-order gravitational theory without the cosmological constant were 
obtained in fllBfl 



Their structure is similar to that of (0), (|32|) and to the solutions (|42|), 
([43|) below. However, it is not possible to pass to a direct limit of conformal invariance in 
the solutions of ||15|| , because the case of conformal invariance is characterized by a nontrivial 
degeneracy, in particular, the massive scalar degree of freedom which is present in a generic 



case is missing here (see also |L4| in this respect). 

Now suppose that a static compact source is composed of identical "atoms" (these may be 
real atoms or elementary particles) and that each of these atoms produces static gravitational 
potential as given by Eq. (|33|) with identical constants m, n, and <fi. In view of the weakness 
of the potential, we also assume the validity of the superposition principle. Then, if fi(r) 
is the spatial distribution of the "atoms" in the source, the total potential is given by the 
expression 



$(r) = J V(\r-r'\)fi(r') d 3 r'. 



(34) 



This potential is the sum of two terms: $(r) = $ m (r) + $„(r). They satisfy the equations 

(35) 



A$ m (r) = 4vrm/i(r), 
A$ n (r) +p$ n (r) = 47rnsin0/i(r), 



that, in the theory under investigation, correspond to the unique Poisson equation of the 
linearized General Relativity. 



For a spherically symmetric compact distribution /i(r), the potential given by Eq. (|34] 
with the kernel given by Eq. ( |33l) is easily calculated: 



M(r') 

r /2 



dr' 



N sin (kr + <fi) Ann sin < 



kr 



■ / /w(r') sin [k (r — r')] r'dr', (36) 

Jr 



where 
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rr 2 Aim f°° 

M(r) = Airni / fi(r')r' dr', N = — — / /i(V) sin (At') r'dr . (37) 
Jo k Jo 

Thus, outside the source, the potential of the form (|33| ) is reproduced with the same phase 
0, but with different coefficients m and n. Moreover, while the coefficient m is additive 
(it plays the role of the gravitational mass of the source), the coefficient n is not: its new 
value N is given by the second expression in Eq. (|37|). However, the coefficient n becomes 
approximately additive for a distribution whose spatial size is significantly less than 1/k. 

If the product kr < 1 in the region of interest (say, on galactic scales), one can expand 
the oscillatory part of Eq. (p3|) in powers of kr to obtain 

V(r) = V - — + — + Qr 2 + O [(kr) 3 } , (38) 

where Vq = — nk cos <j), M = m + nsin0, T = nk 2 sin 0, and Q = q + nk 3 cos 0/6 . We 
thus recover the linear term in the potential of Eq. (|38|), similar to that which was used by 
Mannheim and Kazanas |],[5],|7j to account for the flat galactic rotation curves. However, 
there exists an important observational bound that rules out the possibility for the linear 
term in the expansion (|38D to play a significant role on galactic scales. Note that the 
coefficients — goo( r ) an d g r r( r ) of the metric of our solution are not mutually inverse, which 
is reflected in the fact that the functions a(r) and b(r), given, respectively, by Eqs. (|31~D 
and Q3"2|), are not equal to each other. At small enough distances, both functions reproduce 



Newtonian potentials with the masses, respectively, mo = m + nsin0 and mi = m + 
(nsin0)/2, the difference between them being Am = (nsin0)/2. At the same time, the 
Viking spacecraft observations in the vicinity of the Sun indicate that the ratio Am/m < 2 x 
10~ 3 fl6| (see also 13). This implies the following observational bound for the Sun: 



n sin v . _ _ 3 



rn 



< 4 x 1(T'\ (39) 



Since we assume that the parameter k is sufficiently small so that the expansion (]38[) is 
legitimate on galactic scales, the values of both m and n are additive on such scales and the 
estimate ( j39"D is valid on galactic scales as well. Now, the linear term in Eq. (|38| ) formally 
becomes comparable in magnitude to the Newtonian one only at the distance r ~ \J M/T fa 



m/ (nk 2 sin 0). But, for such distances, we would have kr ~ ym/ (n sin 0) > 10 because of 
estimate (p9|), which contradicts the original assumption kr < 1. Thus, the linear term in 
the expansion (B3) cannot play a significant role on galactic scales, and one should rather try 



the whole potential in the form (^TJ) for a spherically symmetric source with the bound ( |39| ) 
to account for the galactic rotation curves. Such a possibility still remains to be investigated. 

It is instructive to estimate the realistic value of the constant a in Eq. ([!]) for which the 
value of kr is of order unity on a typical galactic scale of 10 kpc, thus making the potential 
of the form (|33|) in principle relevant to the galactic rotation curves. Whatever scalar fields 
are present in the theory, they all contribute to the value of p given by Eq. ([22]). Thus, 
at least the scalar Higgs field of the Standard Model of strong and electroweak interactions 
should be taken into account. The mean value of this field is known to be n ~ 246 Gev, 
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this value will contribute to So in Eq. ( p2]) and, in order that k X (10 kpc) < 1 be valid, we 
must have 

«> 10 74 , (40) 



which, of course, is a severe restriction. It is possible to conceive models in which this 
restriction is weakened; for instance, one can introduce another scalar field with the "wrong" 
overall sign in the action, so that its contribution to the value of the parameter p will be 
of the opposite sign. Its vacuum expectation vaue has to be fine-tuned to counterbalance 
the contribution of rj. Another possibility is to try to construct particle theory without 
fundamental scalar fields, although it may turn out to be very difficult to do this in a 
conformally invariant manner. 

On the other hand, if we take a ~ 1, then the expectation value rj ~ 246 Gev of the 
Standard Model Higgs field leads to the spatial scale 

l/k ~ KT 16 cm (41) 



on which the potential ( j3~3"D oscillates. Its significance might only be manifest on the spatial 
scales of elementary particles, where, of course, the whole theory must be quantized. 

In the case of p < 0, which corresponds to a < 0, the solution for a(r) and b(r) has the 
form 

a(r) = 2qr 2 + m (1 + kr) h n 2 (1 - kr) — , (42) 

b(r) = 2qr 2 + 2n l + 2n 2 — , (43) 



where now k = y/—p , and n\ and n-i are integration constants. Similar solutions in a generic 
second-order gravitational theory (not conformally invariant) without the cosmological con- 
stant were obtained in \TE\. Solutions in the conformally invariant second-order theory with 



the Einstein term but without the cosmological-constant term were also obtained in |14[ 



The physically meaningful solution is selected by imposing boundary conditions at infinity, 
what leads to the condition n 2 = 0. For sufficiently small values of k, the observational 
bound similar to (|39|) implies 

- < 4 x 10- 3 , (44) 
m 



and makes the extra exponential potential in ( pE3"D uninteresting. 

Finally, we note that in the case of p < 0, which corresponds to a < 0, one can also obtain 
solutions by formally replacing the trigonometric functions in Eqs. (|3lD, ([32]), and (|33|) by 
their hyperbolic counterparts and taking k = \f—p- Equations (^) will then remain valid in 
this case as well, with the replacement of sin <p by sinh <p. The structure of the left-hand sides 
of equations (^) reflects, besides the presence of the massless graviton, also the well-known 



S 



presence of a spin-two tachion (in the case of a > 0) or a spin-two massive ghost (in the 
case of a < 0) on the background with S 7^ of the theory described by Eqs. (P, (||) (see 
[p~5| , [T8| j9f1 and references therein). The presence of a tachion in the case of a > indicates 
instability of a large class of classical solutions, including the flat space-time solution in the 
case of A = 0; and the presence of a ghost in the case of a < implies possible absence of 
perturbative unitarity in the corresponding quantum theory. This appears to be the main 
drawback of the conformal theory under discussion. 
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